Preferential affinity of calcium ions to charged phosphatidic-acid surface from a mixed 
calcium/barium solution: X-ray reflectivity and fluorescence studies 
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X-ray reflectivity and fluorescence near total reflection experiments were performed to examine 
the affinities of divalent ions (Ca 2+ and Ba 2+ ) from aqueous solution to a charged phosphatidic- 
acid (PA) surface. A phospholipid (l,2-Dimyristoyl-sn-Glycero-3-Phosphate, DMPA), spread as 
a monolayer at the air/water interface, was used to form and control the charge density at the 
interface. We find that for solutions of the pure salts (i.e., CaCb and BaCl2), the number of 
bound ions per DMPA at the interface is saturated at concentrations that exceed 10 _3 M. For a 1:1 
Ca 2+ /Ba 2+ mixed solutions, we find that the bound Ca 2+ /Ba 2+ ratio at the interface is 4:1. If the 
only property determining charge accumulation near PA were the ionic charges, the concentration 
of mixed Ca 2+ /Ba 2+ at the interface would equal that of the bulk. Our results show a clear specific 
affinity of PA for Ca compared to Ba. We provide some discussion on this issues as well as some 
implications for biological systems. Although our results indicate an excess of counterion charge 
with respect to the surface charge, that is, charge inversion, the analysis of both reflectivity and 
fluorescence do not reveal excess of co-ions (namely, Cl~ or I - ). 

PACS numbers: 



INTRODUCTION 



Varying affinities among same-charge ions to counte- 
rions are expected and are commonly accounted for in 
chemistry. Although they are also common in biological 
processes, they are less pronounced because the strong 
tendency of the aqueous environment in biological sys- 
tems to accentuate the electrostatic interactions rather 
than the chemical ones. However, close to charged inter- 
faces, in channel proteins, or at high salt concentrations, 
the electrostatic description fails at the relevant short dis- 
tances and quantum effects set in. We have undertaken 
the present study to determine the preferential affinity 
of Ca 2+ or Ba 2+ to a highly charged interface formed by 
phosphatidic acid (PA) using X-ray reflectivity and spec- 
troscopy techniques. While PA is a minority lipid in cell 
membranes, it participates in numerous biological pro- 
cesses that are commonly triggered by small changes in 
Ca concentration.— Molecular dynamics simulations 
have predicted that PA has a particular affinity to diva- 
lent ions via its oxygen atoms that may lead to charge in- 
version (namely, the charge density of the bound counte- 
rion exceeds that required to neutralize the surface charge 
exerted by PA) £ These simulations also predict that the 
PA becomes doubly charged by proton release. Recent 
X-ray experiments on PA monolayers have revealed that 
charge accumulation of divalent ions and trivalent ions^2, 
near PA interfaces show charge inversion, i.e. accumula- 
tion of charges that exceed the nominal surface charge. 
It is therefore important for both physical chemistry and 
biology to understand how calcium and other divalent 
ions interact with PA. 
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X-ray scattering techniques to determine 
ion distributions and binding to Langmuir 
monolayers^^^^^^i^i 7 -^^ (LM) have 

been extensively used to address the long standing prob- 
lem of the structure of the electric double layer i 21 i 22 i 23 i 24 
These include X-ray reflectivit y 10 ' 12 and anomalous 
reflectivity , 6 i 7 i 16 i 17 which provide the total number of 
ions and the spatial ionic distributions at the charged 
interface. The near total X-ray reflection fluorescence 
(NTRF) technique monitors specifically and quantita- 
tively ions near the interface ^ 9 ' 12 ' 15 ' 19 ' 20 but does not 
yield spatial information on the counterion layer. The 
near resonance X-ray energy scans at fixed momentum 
transfers can probe specifically the interfacial ions 
and yield their distributions as well as the f'(E) and 
f"(E) fine structures 1 -^ which, like the extended X-ray 
absorption fine structure technique (EXAFS), gives 
information on the local environment of the probed ion. 

While monovalent ion distributions do not show ap- 
preciable specificity, as evidenced by the fact that their 
behavior are well described by PB theory that accounts 
for proton transfer and releas o 21 ' 22 ' 23 from the charged 
headgroup ( 16 i 17 i 18 i 26 distributions for multivalent ions 
show a high degree of ionic specificity.— The theoretical 
understanding of such ionic specificity is not fully un- 
derstood. From classical statistical mechanics, it can be 
shown that Bjerrum pairing can account for some speci- 
ficity related to the different sizes of the ion o 28 ' 29 , but the 
general trends observed in many experiments suggest a 
much stronger specificity, thus requiring more complex 
quantum chemistry calculations. 

Recently Shapovalov and coworkers 20 applied X-ray 
fluorescence studies below the critical angle for total re- 
flection to behenylsulfate monolayers that are spread on 
a variety of mixed ionic solutions to determine preferen- 
tial association of ions to the highly charged sulfate head 
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group. In their study of a 1:1 Ca 2+ /Ba 2+ mixed solution, 
they find that Ba 2+ cations have roughly 10-fold prefer- 
ence over Ca 2+ to associate with the sulfate headgroup. 
Although these authors argue that this result cannot be 
explained by the small difference in the hydration-sphere 
radii, recent theoretical work suggests the opposite 3 ^. 

II. EXPERIMENTAL SETUP AND METHODS 

Ion bulk concentrations were prepared using solu- 
tions of calcium chloride and iodide (CaCl2,CaI 2 ) and 
barium chloride (BaCki) obtained from Sigma- Aldrich. 
Each solution was mixed into ultra-pure water (Milli- 
pore apparatus; resistivity 18.1 Mficm). The surface 
charge density was controlled by a monolayer of 1,2- 
Dimyristoyl-sn-Glycero-3-Phosphate (Monosodium Salt, 
DMPA; MW = 614, from Avanti Lipids, Inc.). Small 
amounts (mili-grams) of DMPA were mixed into a 3:1 
chloroform/methanol solution and spread at the surface 
of the bulk concentrations in a temperature controlled 
Langmuir trough. The trough is placed in a sealed can- 
ister that is kept at a constant temperature (22°C) and 
is continuously purged with water-saturated helium gas. 
10-15 minutes were allotted for solvent evaporation af- 
ter which the monolayer was compressed at a rate of 1 
A 2 per molecule per minute to the desired pressure (the 
experiments reported here were conducted at a constant 
pressure ir rs 25mN/m). During the monolayer compres- 
sion the surface pressure was recorded by a micro-balance 
using a filter-paper Wilhelmy plate.— 

X-ray reflectivity and fluorescence studies of monolay- 
ers at gas/water interfaces were conducted on a home 
built liquid surface X-ray reflectometer, using an Ultra- 
X18 Rigaku X-ray source generator with a copper ro- 
tating anode operating at ~ 15 kW power. Down stream 
from the beam source, the characteristic Cu Kq (A = 1.54 
A) is selected and deflected onto the liquid surface at a 
desired angle by a Ge(lll) crystal. Scattered photon in- 
tensities were normalized to an incident beam monitor 
placed in front of the sample. The X-ray diffractometer 
was used in two configurations: reflectivity and fluores- 
cence. 

Specular X-ray reflectivity experiments yield the elec- 
tron density (ED) profile across the interface, and can 
be related to molecular arrangements in the film . 31 ' 32 ' 33 
The reflectivity for a LM at Q z is calculated by 

R(Qz) = Ro(Qz)e- Q ^ , (1) 

where Rq(Q z ) is the reflectivity due to an ED profile 
composed of step functions, calculated by the recursive 
dynamical method,— and a is an effective surface rough- 
ness, which has a weak Q z dependence according to capil- 
lary wave theory applied to X-ray reflectivity studies ! 35 ' 36 
In our reflectivity measurements, we use variable slit con- 
ditions that roughly keep the incident beam footprint 
constant on the liquid surface and also adjust the detec- 
tor slit accordingly. This increases the angular accep- 



tance of the detector with Q z , which makes the depen- 
dence of a on Q z large enough to require a correction. 
To a first order approximation, we correct for that by 
assuming a 2 = a + bQ z , where a and b are two fitting pa- 
rameters. Other variable parameters used to construct 
the ED across the interface include the thicknesses of the 
slabs (di) and their corresponding electron densities (pj) 
(for details on the definitions of parameters see Ref . \17l ) . 
The minimum number of slabs is the one for which the 
addition of another slab does not improve the quality of 
the fit, i.e. x 2 > significantly. 

To collect the fluorescence signal from the monolayer- 
counterion system, we use an energy dispersive detec- 
tor (X-PIPS Detector; model SXP8, with a Multiport 
II multichannel analyzer both from Canberra Industries 
Inc.). The pencil- like tip of the detector is lowered to the 
surface in an aluminum well in front of a thin Kapton 
window located about ~ 2 cm above the liquid surface 
for fluorescence measurements (the fluorescence-setup is 
under the same conditions as for reflectivity, namely, the 
canister holding the trough is sealed and is purged with 
water saturated helium). The major advantage of the 
fluorescence technique is that it can distinguish between 
different elements by their characteristic emission lines. 
In the present study due to the low-yielded signal, the 
fluorescence signal is distinguished and displayed over 
two regions corresponding to two regions of the incident 
beam angles: below the critical angle for total reflection 
(Q c ~ 0.022 A -1 with little variation due to salt in the 
solutions), which is highly sensitive to elements at the 
surface, and above the critical angle, which is sensitive 
to elements in the bulk and surface of the solution. Below 
the critical angle, the evanescent X-ray waves penetrate 
to 60 — 80 A from the surface, which makes the bulk con- 
tribution to the fluorescence intensity almost negligible 
at concentrations smaller than 10~ 3 M. 



III. EXPERIMENTAL RESULTS 



A. Isotherms 

Figure Q] shows surface pressure versus molecular area 
isotherms for three different salt solutions (as indicated) 
at concentrations of 10 _3 M with DMPA and pure wa- 
ter with DMPA. The isotherm of DMPA on pure water 
lacks any distinct phase transitions with a smooth in- 
crease with compression. For DMPA spread on salt so- 
lutions, two slopes associated with crystalline tilted and 
untilted phases of the hydrocarbon chains with respect 
to the surface normal can be clearly identified. As the 
monolayer is compressed, it reaches the untilted stage 
(tt ~ 13 mN/m) at the sufficiently low molecular area 
of - 40 A 2 per molecule. The isotherms of DMPA on 
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Figure 1: Surface-pressure versus molecular area (n — A) 
isotherms of DMPA spread on water and dilute solutions, 
10~ 3 M CaCl 2 , BaCl 2 , and a mixture solution (5 x 10" 4 M 
CaCl 2 + 5 x 10" 4 M BaCl 2 ). Although the isotherms on salts 
solution are different than that of water they are hardly dis- 
tinguishable with no specific features that can be associated 
with the different ions. 



the different salt solutions vary just slightly, exhibiting 
little or no specific features to the different ions in the 
solutions. All fluorescence and reflectivity measurements 
were conducted in the untilted crystalline stage (it ~ 25 
mN/m). 



B. X-ray Reflectivity 

Normalized reflectivity data (R/Rf, where Rf is the 
calculated reflectivity of pure water with an ideally flat 
interface; a = 0) of DMPA monolayers spread on 1CP 3 M 
solutions and on pure water are shown in Fig. [2j It is 
evident from Fig. that the overall reflectivities for the 
samples prepared on salt solutions are significantly higher 
than that of the pure water data. This is a qualitative 
evidence that the ions (Ba 2+ and/or Ca 2+ ) are accumu- 
lating at the surface. Also, there is a visible shift in the 
minima of the reflectivities to a smaller Q z for the solu- 
tions compared to water reflectivity, indicating a thicker 
layer due to the association of the ions at the headgroup 
region. 

The normalized reflectivity data were fit (solid lines 
Fig. [2ji and b) to a two slab model, one associated with 
the hydrocarbon chains of DMPA close to the vapor, 
and the other, contiguous to the bulk solution, associ- 
ated with the phosphatidic-acid headgroup and bound 
ions and water molecules. The fitting parameters for the 
reflectivity data are listed in Table HI Using these pa- 
rameters the electron density profiles for each solution 
and for pure water are obtained and plotted in Fig. [2fc. 
It is interesting to note that the electron densities for 
divalent ion sub-phase are, within experimental uncer- 
tainties, identical. Qualitatively, this surprising result 
indicates there must be more calcium ions per DMPA 
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Figure 2: (a) Normalized reflectivity data (R/Rf) versus mo- 
mentum transfer Q z for each of the dilute solutions (as indi- 
cated) and on water with DMPA on a linear scale, (b) The 
same data on a logarithmic scale (curves and data are shifted 
by a decade each for clarity), (c) Electron density profiles 
used to calculate the best fits shown by solid lines in (a) and 
(b). The electron densities are generated with the parameters 
listed in Table H 



than barium ions as the two ions, Ba 2+ and Ca 2+ , have 
very different numbers of electrons, 18 and 54, respec- 
tively. The phosphatidic-acid headgroup on pure water 
has a smaller thickness and smaller electron density than 
those obtained for the salt solutions, as expected since it 
does not have an electron-rich ion layer as with the salt 
solutions. 

From the structural parameters in Table HI one can de- 
rive the number of ions per DMPA molecule at the surface 
by applying volume constraints and assuming the co-ions 
(e.g., CI - ) do not penetrate into the headgroup region (as 
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Table I: Parameters used to fit the reflectivity data as well 
as the derived number of ions and water molecules at the 
interface. Uncertainties are obtained by fixing a parameter at 
values away from the optimum and readjusting all of the other 
parameters to a new minimum until \ 2 increases by 25%, and 
by imposing limiting values on some of the parameters to keep 
them within physical range; for instance, the electron density 
of the hydrocarbon chains was kept within 0.32 ±0.03. 

shown below). We use the following equations to solve 
for the number of water molecules (n w ) and the num- 
ber of ions (n,) per phosphatidic acid in the headgroup 
slab, 3 ! 

Ad h = V h + n w V w + mVi (2) 

Ad h ph = N h + n w N w + n t N % , (3) 

where A is the molecular area obtained from the isotherm 
and diffraction data,— i 31 i 37 dh and ph are the thickness 
and electron density of the head group, respectively, 
Ni is the number of electrons per molecule (for wa- 
ter molecule N w = 10, and for Ca 2+ and Ba 2+ 18 
and 54 electrons, respectively). The molecular area 
can be confirmed self consistently by the electron den- 
sity of the hydrocarbon chains p c hain from the relation 
A = N cha i n /{d chain Pchain)- The volume of the phospha- 
tidic headgroup, Vh, is obtained by the relations above, 
eqs. [2] and [3j using the the electron density profile pa- 
rameters for DMPA spread on pure water. The precise 
position that parses the molecule into the hydrocarbon 
chain compartment and the headgroup region influences 
just slightly the derived parameters, and this change is 
accounted for in the uncertainties of the parameters given 
in Table [1 As seen in the table and as expected by the 
electron-density profiles, the number of Ca 2+ and the 
number of Ba 2+ ions per DMPA is not the same but 
rather seems to favor Ca 2+ by a ~ 2 : 1 ratio for the 
unmixed salt solutions. 



C. Fluorescence 

Fluorescence scans were performed on all the solutions 
with and without DMPA, and the background data was 
obtained from a pure water sample without DMPA and 
subtracted from all data to remove electronic and stray 
signals from the trough. The data from the fluorescence 
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Figure 3: (a) Fluorescence signals below the critical angle 
for 10~ 3 M BaCb of the solution only and the same solution 
with the DMPA monolayer as indicated, (b) Fluorescence 
data above the critical angle from the same two samples. The 
signals above and below the critical angle were integrated over 
a Q z range as indicated in the figures. 

were then integrated over Q z ranges to improve signal 
to noise ratio as described above. Figure [3] shows fluo- 
rescence data below (a) and above (b) the critical angle 
for 1CT 3 M BaCl 2 with and without DMPA, respectively. 
The intensity was integrated over the range shown in the 
figures. As shown in Fig. [3](a) (below the critical angle), 
without a monolayer the evanescent X-ray beam does not 
penetrate far enough into the bulk to generate any visible 
fluorescence intensity. However, with the DMPA mono- 
layer on, the ha (4.46 keV) and L/3i (4.8 keV) emission 
lines for Ba 2+ are clearly visible, which indicates that 
the concentration of ions at the interface is much larger 
compared to that of the bulk concentration. In Fig.[3](b), 
emission lines can be seen in both the solution without 
and with DMPA. However, the solution with the mono- 
layer has an enhancement for each emission line as the 
ions at the surface contribute significantly to the overall 
intensity. The difference in the intensities of the two data 
sets in Fig. [31(b) quantifies the contribution of the ions at 
the surface. Spectra of Ka « 3.7 keV and Kf3 « 4.03 
keV emission lines for the CaCl2 solution (not shown) 
essentially show similar trends. Figure [4] shows that 
the emission lines below the critical angle (contributed 
primarily by surface) are practically identical at 10 _3 M 
and 10 _2 M BaCl 2 , indicating that the monolayer is satu- 
rated with ions for concentrations stronger than 10 _3 M. 
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Figure 4: Integrated fluorescence signals below the critical 
angle for solutions of 10~ 2 M and 10 -3 M BaCb showing the 
binding of ions reaches saturation for concentrations < 1CP 3 
M. 



CaCl2 and the mixture solutions display similar binding 
phenomena at these concentrations. 



To get more quantitative results we define a "standard 
ratio" as the ratio of fluorescence signals for specific emis- 
sion lines from two different elements in equal quantities 
in the bulk solution (e.g., Ca 2+ /Ba 2+ ). Prom Fig. Ufa), 
pure solutions without DMPA, we establish a "standard 
ratio" (by comparing the peak intensities of Ba 2+ La and 
Ca 2+ Ka) of the emission lines generated by equal con- 
centrations of Ba 2+ and Ca 2+ ions at ~ 6 : 1 in favor 
of Ba 2+ . Comparing the ratio between the same peak 
intensities provides the relative population of each ion at 
the interface. From Fig. EJb) we find the Ba 2+ /Ca 2+ 
ratio of the pure solutions reduces to ~ 3 : 1 and for the 
mixed solution the ratio is even lower ~ 1.5 : 1. This 
shows that Ca 2+ in unmixed solution populates the PA 
interface twice as much as the Ba 2+ does, consistent with 
X-ray reflectivity results. In the Ba/Ca 1:1 mixed solu- 
tions calcium ions are more likely to associate with the 
PA at the interface by a ratio of ~ 4 : 1 over Ba 2+ . 



IV. DISCUSSION AND CONCLUSIONS 

The objective of this study has been to explore the 
binding specificity of divalent ions in solutions to a phos- 
phatide acid charged surface and in particular, to deter- 
mine which of the two Ba 2+ or Ca 2+ is more favorable to 
binding to the DMPA monolayer. Our results show that, 
in a 1:1 mixed solution, Ca 2+ is favored over Ba 2+ by 
almost a 4:1 ratio. This is strikingly different from the 
results of Shapovalov and coworkers^, who found that 
Ba 2+ outnumbers Ca 2+ by a 10:1 in binding to a be- 
henylsulfate monolayer. Clearly, the only possible origin 
of such a discrepancy is a strong specificity of ionic and 
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Figure 5: (a) Integrated fluorescence signals above the critical 
angle showing the relative intensities of Ba and Ca emission 
lines from the solutions only at 10~ 2 M (no monolayer). This 
data is used to calibrate the relative intensities from the same 
number of Ba and Ca ions, (b) Integrated fluorescence sig- 
nals below the critical angle for DMPA monolayers spread on 
BaCb, CaCi2 and 1:1 mixture at a nominal 10~ 2 M. This 
indicates the amount of bound Ba 2+ to PA is significantly 
smaller than that of Ca 2+ for DMPA on the mixture. 



headgroup interactions, and cannot be quantitatively ac- 
counted for by theories based on classical electrostatic in- 
teractions, even if the ionic size is properly accounted for. 
Furthermore, our finding that Ca 2+ binds more strongly 
than Ba 2+ can only be understood if the ion dehydrates 
upon binding. Quantitative theoretical account for these 
effects remains an outstanding theoretical problem. 

Our results also imply charge inversion at the surface in 
particular, for the CaCl2 and Cal2 solutions, as shown be- 
low. This happens when the charge density of the coun- 
terions exceeds the surface charge density exerted by the 
phosphatidic monolayer. Our analysis shows there there 
is one Ba 2+ per DMPA molecule, which is just enough 
to neutralize the surface charges (doubly charged phos- 
phatidic acid), whereas there is more than one Ca 2+ ion 
per DMPA molecule at the surface. This suggests charge 
inversion, from negatively charged to positively charged 
interface. In the fluorescence experiments above, we did 
not detect strong emission signals from Cl~ to suggest 
that there exists an accumulation of the co-ion (Cl _ ) at 
the interface due to charge inversion. To further explore 
if the anions are rising to the charge-inverted surface, we 
performed another experiment using 10 -2 M Cal2 solu- 
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Figure 6: (a) Fluorescence signals from the surface of Cal2 
(10~ 2 M) solution only, and the same solution with the DMPA 
monolayer on. At this relatively high concentration the signal 
from iodine is clearly seen before spreading the monolayer 
whereas the signal form calcium is very weak to detect. With 
the DMPA monolayer on, the Ca signal is enhanced whereas 
the signal from the iodine slightly weakens. This shows there 
is no enrichment of the iodine co-ion at the interface due to 
possible charge-inversion, (b) The electron density profiles for 
films of DMPA spread on 1(T 2 M Cal 2 and 10" 2 M CaCl 2 
solutions. This shows the amount of I~ is negligible at the 
interface although the number of Ca 2+ ions exceeds those 
necessary to neutralize the phosphatidic-acid. 



tion with DMPA as a monolayer to improve the fluores- 
cent signal from the co-ion, in this case from I - . Figure 
EJa) shows spectra below the critical angle, indicating 
surface enrichment of Ca 2+ but no significant difference 
in the emission lines of I - before and after the spreading 
of DMPA. We have also measured the reflectivity from 
the same DMPA monolayer on Cal 2 solution to detect 
any changes in the electron density near the headgroup 
region compared to that of DMPA on CaCl2 solution, 
that may indicate the accumulation of I~ at the inter- 
face. Figure E](b) shows the ED profiles of the films on 
both solutions (CaCl2 or Cal^) are practically the same. 
We therefore conclude that if there are co-ions (iodine or 
chlorine ions) at the surface due to charge inversion, their 
concentration is below 0.1 iodine ion per PA molecule ac- 
cording to our detection sensitivity, or that the Ca 2+ ions 
migrate to the surface partially as CaOH + . 

To summarize, in numerous physiological studies such 
as those of the 'Ca receptor' in preganglionic nerve ter- 
minals (interacting primarily with acetyl-choline) it has 
found that Ba 2+ ions, do not appear to compete with 
Ca 2+ for the receptor—. This is consistent with our stud- 
ies that show stronger affinity of Ca over Ba to negatively 
charged PA interface. Systematic studies of model sys- 
tems using spectroscopic X-ray techniques that are ion 
specific open the door to unraveling the mechanism of 
this and related physiological phenomena. 
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